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ie I. INTRODOCTION 


A. GENERAL DESCRIPTION 


Since the implementation of the Finite Element Analysis 
Program (hereafter refered to as FEAP) at the Naval 
Postgraduate School, it has had only limited use. It is now 
available to general NPS users in an interactive form 
through the IBM 3033 VM/MVS time sharing systen. The data 
management and user interactive command structures were 
established within the existing interactive progran, 
Computer Analysis Language (CAL). This integration of 
systems provides the ability to solve linear and nonlinear, 
steady and unsteady, two and three jimensional heat conduc- 
tion problems involving temperature dependent thermophysical 
properties and complicated radiation/convection bourdary 
conditions. 

Additional capability was provided +o CAL system users 
both in changes to existing subroutine groups and in the 
addition of a new graphics group with its attendant command 
structure. The graphics functions enable the user to plot 
two- and three- dimensional structural and heat transfer 
mes hes. Another new facility is the HELP operation which 
allows a user experiencing trouble with a particular opera- 


tion to interactively obtain assistance. 


Be. HISTORICAL BACKGROUND 


The original CAL program was developed by Professor E. 
Le. Wilson of the University of California in 1977 [Ref. 1}. 
It was later adapted and modified for use at the Naval 
Postgraduate School by LCDR L. B. Elliot [Ref. 2]. The FEAP 
program -wasS written by Professor R. ite Taylor of the 
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University of California [Ref. 3] in 1977. Implementation 


of FEAP at Naval Postgraduate School was done by LT J. M&M. 
Bettencourt [Ref. 4]. 


C. OBJECTIVES 


The objectives of the author's work have been to: 


1. 


integrate the data management system of FEAP with CAL 
to create an interactive conduction heat transfer 
problem solving system; 

modify existing operations to extend their useful- 
ness; 

add operations to extend the capabilities of the 
program to include graphics; 

create a HELP facility; 

create a USER'S MANUAL to facilitate use of this 
progran. 
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II. ORGANIZATION OF HEAT TRANSFER GROUP 





This chapter provides a general overview of the organi- 
zation of the Heat Transfer Group of CAL. It is intended to 
provide sufficient information to permit users to operate 
the CAL heat transfer package. 

The execution of the program is flexible and controlled 
by user selection of operations in a logical sequence from 
the commands that are available. There are two broad 
categories of operations, data input and problem solution. 


A. DATA INPUT 
1. Initialization 


The heat transfer yroup performs matrix creation and 


Manipulation automatically. As the problem progresses, 
arrays are created, altered and deleted under progran 
control. Through the HTXFR operation the user provides 


sufficient information to establish the initial arrays for 
data input and problem solution. The number of nodes, 
number of elements, number of material sets, spatial dimen- 
sion, number of degrees of freedom per node and the maximun 
humber of nodes per element are required. For heat transfer 
problems the number of degrees of freedom per node is always 
one. The option to assigna higher number is available 
because the equation solvers in this program are applicable 
to other fields of which future work may make use. 


2. Nodal Coordinates 


Nodal coordinates are input via the COORD operation. 
This operation has built in node generation capability. BY 


specifying an initial point and a node generation vector, 
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the user may easily input large meshes. Coordinate system 
conversion is also available. Coordinates may be input in 
the Cartesian system, the cylindrical system with any one of 
the three axes longitudinal, the spherical system or any 
combination of the above systems. All coordinatesS are 
converted to Cartesian coordinates for use with CAL. 


The ELCON operation inputs the element connectivity 
data. Here again is a generation capability. The user may 
specify the connectivity for one element and a generation 
vector to create additional rows. 

For the two-dimensional elements the user may 
specify a 4 to 9 node isoparametric element. There is an 8 
to 21 node isoparametric element for three-dimensional 
elements. Both of these elements must follow the numbering 
convention shown in Appendix A. 


GY. Material Properties 


The required amount of material property information 
varies from problem to problen. The PROP operation prompts 
the user for the information required to solve the problem 
at hand. At a minimum the material's conductivity (k), 
specific heat (c), specific mass (rho),heat generation per 
unit volume (q''**) and the geometry type (plane or axisyme- 
tric) must be specified in a consistent system of units. 
Appendix A includes examples of consistent systems. 

The user also inputs the number of GauSSian points 
per direction for quadrature and codes for temperature 
dependent properties and boundary conditions. The codes 
indicate other information required. Temperature dependent 
properties are input as tables and linear interpolation is 
used to determine the property value at a given temperature. 
Boundary conditions are specified as shown in Appendix A. 
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If an exterior boundary line or surface condition is not 


specified, it is assumed to be insulated. 
5- Constant 


For problems requiring certain nodes ‘to be at 
constant temperature, the CTEMP operation is available. 
This operation may also be used during the solution stage of 
the problem to provide step changes at previously specified 
constant temperature nodes. Because this operation gener- 
ates an array used in profiling the solution equations, the 
user may not change the node numbers that were established 
as constant constant temperature nodes after execution of 
the PROF operation. The temperatures of these nodes, 
however, may be changed. 


6. Equation Profile 





The PROF operation establishes the equation profile 
for problem solution. Prior to the execution of this opera- 
tion, any data input may be changed by specifying the appro- 
priate operation and re-entering the data. After its execu- 
tion the user may not change the nodes designated as 
constant temperature nodes to temperature varying nodes or 
vice versa. The user may change the value of the constant 
tem peratures. 


Be SOLUTION 


The matrix formulation of the heat transfer problem as 


discussed in Reference 4 is: 
(K]{T} + (C}{T} + {(F} = (0} (1) 


where (K) represents the conductivity matrix, (C) represents 


the heat capacity matrix, {F} represents the flux vector and 
{T} represents the temperature vector. The derivative of 


{T} with respect to time is {fT}. The flux vector includes 
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heat generated per unit volume and boundary fluxes as speci- 
fied for given boundary surfaces. This is a fully general- 
ized formulation, including non-linearities, since the 
matrices {K] and [CJ] and the vector {F} can be temperature 
dependent. 


For heat transfer problems, which involve only one 
degree of freedom per node, the conductivity matrix, [K], 
Will always be symmetric. The capability for generating 
uns ymmetric matrices was provided, but will only be appli- 
cable when additional types of problems are programmed into 
the CAL systen. The command for the unsymmetric conduc- 
tivity matrix formulation is USYMC. 


2. Forming Heat Capacity Matrix 


The heat capacity matrix, [CJ], used in time depen- 
dent problems, can be generated with either CCAP or LCAP. 
To form a consistent capacitance approximation use the oper- 
ation CCAP. A lumped capacitance approximation is formed 
uSing the LCAP operation. 


3- Forming Flux Vector 


To complete the problem formulation, the flux 
vector, {F}, must be generated. The FORM operation forms 
the flux vector taking into account the internal heat gener- 
ation and boundary surface fluxes as indicated in the PROP 
operation. 


4. Equation Solving 


Once the time independent problem is formulated, the 
temperature vector, {T}, is calculated by the CALC opera- 
tion. Time dependent problems do not use this operation, 


but rather the ordinary differential equation solver. 
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Ordinary Differential Equation Solver 


The first order ordinary differential equation 
solver is accessed with the ODE operator. It employs the 
Zienkiewicz two- and three-level schemes [Ref. 4, 5]. 

In addition to the time step size change using the 
DTIM operation, an optional automatic time step adjustment 
is incorporated in the ODE operation. The norm of the 
difference between temperature vectors at two consecutive 
times is computed at each step. If the norm is less than a 
user specified maximum temperature difference, the time step 
is doubled before going to the next step. If the norm is 
greater than a user supplied minimum temperature difference, 
the time step is halved and calculation for that time step 
is repeated until the norm is acceptable. If the tempera- 
ture differences are specified as zero, no time step adjust- 
ment will be performed. 

The user specifies one of three functions which are 
performed by the ODE operation. A second operation name, 
which must be separated by a comma, follows ODE. The 
options are INIT, LINE or QUAD. 

The operation ODE,INIT is used to input the integra- 
tion constants theta, beta and gamma [Ref. 5], the maxinun 
and minimum temperature differences for the automatic time 
step adjustment and the initial temperature vector. No time 
integration is performed by this instruction. 

ODE,LINE performs the two point scheme and the 
current temperature vector is substituted by the newly 
calculated temperature vector. 

The ODE,QUAD operation is similar to the ODE,LINE 
Operation but uses the three point scheme. 


6. Printing Nodal Temperatures 


Once the temperature vector has been updated, the 
PTEMP operation prints the temperatures in node number 
order. 
16 





III. ORGANIZATION OF GRAPHICS GROUP 


The graphics group is capable of displaying meshes 
either on the IBM 3277 dual screen terminal system or any 
PLOT-10 compatible terminal. It is initiated through the 
use of the GRAPH operation by which the user specifies the 
type of graphics terminal in use. 


A. TITLE 


The TITLE operation is used t95 title the mesh being 
displayed. The user may input up to three lines of fifteen 
characters for the title. This operation must be specified 
immediately prior to the displaying operation. 

This operation calls the USRIN subroutine which reads 
from the terminal three lines of characters. A flag (IFLAG) 
is set to indicate a title is to be printed. 


Be HEAT TRANSFER MESH 


The PLHX operation locates the heat transfer coordinate 
and connectivity arrays and displays the mesh. The viewing 
area is optimized so the longest dimension is full screen. 
The maximum and minimum coordinate values for each direction 
are displayed. This operation displays two- and three- 
dimensional meshes. If the mesh is three-dimensional, «he 
user must specify the viewing plane. 

This operation calls either the FPPLOT (2-D) or FP3PLT 
(3-D) subroutine. These subroutines scan the coordinate? 
arrays for maximum and minimum values and initialize the 
graphics screens. They both use subroutine BOX to set the 
virtual window, set the screen window and draw a bex around 


the plotting area. The standard element connectivity is 
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stored in an array which is used in conjunction with the 
user's element connectivity matrix to draw the elements line 
by line. The center node is plotted using a “+#" symbol. 
The user's title is plotted (if IFLAG is not 0) by subrou- 
tine USRTIT, which then resets IFLAG. The maximum and 
Minimum values of the vertical and horizontal coordinates 
are printed by SCRDAT. Subroutine SCRDAT also indicates axis 
orientation. Prior to terminating the screen graphics, 
subroutine TITLE is used to write the title box identifying 


the type of mesh being displayed. 


C. STRUCTURAL MESH 


The PLST operation displays the structural system speci- 
fied by user supplied coordinate and connectivity arrays. 
The viewing area is optimized as described for the PLHX 
operation. Likewise, the user may specify the viewing plane 
as either the X-Y, Y-Z or X-2Z2 plane. 

This operation calls the CLPLOT subroutine. CLPLOT uses 
the same logic and subroutines as FPPLOT. The difference 
between the two subroutines is in the method of storage of 


the coordinate arrays and the connectivity matrix. 


18 





IV. CHANGES TO EXISTING GROUPS 


Modifications were made in several existing subroutine 
groups to make their operations more versatile. One utility 
subroutine was improved and one was added. 

The changes that were made are sensitive to prior 
versions of CAL. The same results are obtained for previ- 
ously existing operation commands. No files used with other 
editions of CAL need to be modified to operate with this 


progran. 


A. UTILITY SUBROUTINES 


1, Subroutine RCAR 


Subroutine RCARD reads and interprets the operation 

commands. An operation command has the form 
OP,M1,M2,M3 ,44,M5 ,N1,N2,N3,N4 

where OP is the operation name, M1 to M5 are matrix names 

and N1 to N4& are integers. Previously, N1 to N& had to he 

values greater than or equal to zero. The symbol/state 

logic matrix and subsequent action codes were modified to 

allow users to input negative integer values. 

This is important. for the ZERO operation which 
enables users to create matrices with a given value in the 
diagonal locations and another value in the off diagonal 
locations. Prior to this change, a user desiring to create 
such a matrix with negative values had to input the matrix 
row by row. 


2- Subroutine FRTCMX 


Subroutine FRTICMX is a new utility subroutine which 


allows the CAL program to invoke most CP/CMS commands. 





After the invoked command is executed, control returns to 
CAL. 
It is presently used in the SAVE and RESUME opera- 

tions to invoke the CMS command 

FILEDEF NSAVE DISK M1 SAVE (RECFM VS LRECL 7290 

BLKSIZE 7294) 
where NSAVE is a logical unit number assigned to the SAVE 
and RESUME operations and M1 is a user input file name. 


Be. GROUP 1 


1. SAVE Operation 





The save operation creates a file on the user's 
A-disk containing all arrays in storage at the time of 
issuance. Previously the entire 100,000 word main array was 
stored, regardless of how many words were actually being 
used. It was always stored under the name FILE 02, 
preventing the user from saving more than one problen. 

The method of storing the array was changed to store 
just the number of locations actually being used. More than 
one problem may be saved because the subroutine FRTICMX was 
used to create a SAVE file with a name assigned by the user. 
If a name is not specified FILE 02 will be the name of the 
saved problen. 
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mo 


ESUME Operation 


The RESUME operation reads a saved file into memory. 
It was altered to read named files saved by the new SAVE 
Operation. If a name is not specified FILE 02 will be read. 


3. 


IS 
I" 


A Operatio 





The LOADI operation loads integer arrays. The arays 
were input row by row. The option to generate arrays was 


added. The user may specify one row and a row generation 
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vector. The number of rows specified will be automatically 
generated. This operation was moved from the static anal- 


ysis subroutine group to the general matrix command group. 


4. PRINT Operatio 





The PRINT operation prints an array in matrix 
format. Previously it could only print arrays containing 
real numbers, resulting in the erroneous printing of arrays 
containing integers. The user may now specifiy whether the 


array to be printed contains real or integer values. 


5. HELP Operation 


The HELP operation was added to the general command group 
and provides the user with information on the use of all the 
available operations. It accesses a file of instructions, 
sorts through them to find the desired operation and 


displays the appropriate information on the screen. 


C. GROUP 2 


1. NODE 


Operation 





The NODES operation creates the matrix of nodal 
coordinates for a structural problen. This information was 
entered in cartesian coordinates, node by node. 

The user may now opt to enter data in cartesian, 
cylindrical (any axis longitudinal) or spherical coordi- 
nates. It will be converted to the cartesian coordinate 
System used by CAL. Additionally, the user may generate new 
nodes by specifying one node and a node generation vector. 
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V. SOLUTION OF CONDUCTION HEAT TRANSFER PROBLEMS 


In order to solve a conduction heat transfer problen, 
the user must provide the solution algorithm to CAL. A 
discussion of the matrix manipulation and equation solving 
techniques can be found in Reference 4 . This chapter will 


present possible algorithms for solution of this class of 
problems. 


A. STEADY STATE PROBLEMS 
These problems take the form: 


CK){T} + (F} = {0} (2) 


1. Linear Heat Conduction Problems 


This is the simplest case to consider. After input- 
ting mesh data, the user must form the conductivity matrix 
(K], (SYMC) and the load vector {F}, (FORM). Then the nodal 
temperatures must be calculated (CALC) and printed (PTEMP). 

Consequently the segquence of solution operations for 
this type of problem would be: 
oan C 
FOR 
CALC 
PTEMP 


2- Non-linear He 


c+ 
Oo 
12 
Qu 
S 
(A 
lh 
O 


Since the conductivity matrix is time dependent on 
temperature, an iterative algorithm must be used. This 


requires a looping operation (LOOP, NEXT) around the linear 
Steady state sequence. 


me 





The solution operations may be: 
LOOP, N1 
SYMC 
FOR M 
CALC 
PTEMP 
NEXT 
where N1 is the user's guess of the number of iterations 
necessary to obtain equilibrium. However, the program main- 
tains an internal check on the residuals. When they 
decrease below the predefined tolerance (TOL command, 
default is 10-9), the looping operation is *erminated upon 
the subsequent NEXT command. 


Bo. TIME DEPENDENT PROBLEMS 


These problems involve the full form of equation (1): 


CK}(T} + (C}{T} + (F} = 0} (1) 


1. Linear Heat Conduction Problems 


This case requires the solution of a first order 
ordinary differential equation (ODE). Additionally a heat 
capacity matrix, [C], must be formed (CCAP or LCAP) anda 
time step provided (DTIM). 

The differential 2quation solver is accessed using 
the operation ODE,M1 where M1 is one of the following: 

INIT to specify initial temperature vector and 
the integration constants 

LINE +o perform the two-time level algorithn 

QUAD +o perform the three-time level algorithn 

The heat capacity and conductivity matrices are 
unchanged in a linear problem and must be placed outside the 
loop. The load vector is reformed every time step and the 
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time must be advanced with ADTIM. These operations are 
included in the loop. 
The sequence of operations to solve a linear time 
dependent problem may be: 
DTI M 
ODE, INIT 
SYMC 
CCAP (or LCAP) 
LOOP, N1 
FOR M 
ODE,LINE (or QUAD) 
ADT IM 
PTE MP 
NEXT 
where N1 is the number of time steps th2 user wants to take. 


2. Non-linear Heat Conduction Problems 


The heat capacity and/or conductivity matrices are 
temperature dependent in this class of problem. The temper- 
ature dependent matrix (matrices) must appear inside the 
looping operation whereas the constant property matrix would 
be excluded from the loop. 

To solve a fully non-linear problem the following 
Operation sequence may be used: 

DTI M 

ODE, INIT 

LOOP, N1 

Sy c 

CCAP (or LCAP) 

FOR M 

ODE,LINE (or QUAD) 
ADT IM 

PTE MP 

NEXT 


where N1 1s the number of time steps the user wants to take. 
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C. NOMERICAL EXAMPLES 


A hollow cylinder with a four inch outer diameter 
and a three inch inner diameter was subjected to an axial 
forced connvection condition in awind tunnel by Professor 
P. F. Pucci of the Naval Postgraduate School. There were 
sixteen one-quarter inch wide heating strips equally spaced 
over 180 degrees of the outer surface as illustrated in 
Figure 1. 

The heating strips were maintained ata constant 
temperature of 1609 F and measurements of the surface temp- 
eratures between the strips were made using teledeltos 
paper. The ambient temperature was 60° F. The cylinder was 
considered to have a density of 70 lb/ft? and a specific 
heat of 0.6 BTU/1lbm°F. Tests were made with four heat 
transfer coefficients (h). 

The finite element model tock advantage of the 
cylinder's symmetry, consisting of one half of a heating 
strip plus one half of the interval between strips. The 
element mesh,as generated by the PLHX operation, is shown in 
Figure 2. The comparison of the model data to the tele- 
deltos paper measurements was very favorable and is shown in 


Figure 3. 


2. Transient Surface Temperatures in an Infinite Plate 


EEE = WES ea «> ae ‘=a 


An infinite flat plate was considered as a test 
problem (Figure 4). The plate was at an initial uniform 
temperature greater than the ambient temperature and then 
exposed to convection conditions with a constant external 
heat transfer coefficient (h). 

The approximate solution was obtained using a 
Heisler chart. The temperature of the outer wall was 
computed at 1 minute, 5 minutes and every 5 minutes there 


25 





after up to 30 minutes. The comparison between these values 
and those generated by CAL was very close and 1S shown in 
Figure 5.~ 


26 


reelavszmir weree?d noes 
Te eb" : 2 ePenie "Tae — 











VI. CONCLUSIONS AND RECOMMENDATIONS 


The code that was integrated into CAL provides an accu- 
rate and reliable means for solving a variety of conduction 
heat transfer problems. The system is user friendly both in 
prompting for input and detecting errors. The use of the 
heat transfer group of commands is encouraged, as well as 
efforts to increase its versatility. 

While the capabilities of the program are significant, 
there is room for improvement. The present version uses 
primarily an in-core solution technique, which restricts the 
problem size to within the user's virtual machine space. 
The capacity to handle larger problems may be increased 
through the use of an out of core technique for building and 
storing the conductance and capaitance matrices, as well as 
equation solving. . 

The graphics package could be expanded and improved to 
provide more information to the user. Ther plotting 9 Of 
isoparametric elements needs to include the capability to 
generate and plot curved lines. Another desirable capa- 
bility is to portray all three axes on the screen with 
apparent depth and allow rotation to any desired view. A 
capabilty to plot isotherms could be added to enhance the 
interpretation of the results. 
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APPENDIX A 


USER'S GUIDE 


This appendix provides details on the use of CAL with 
the IBM 3033 computer at NPS. The program as modified at 
the Naval Postgraduate School is subsequently referred to as 
CAL-NPS. Section A provides the details on the command 
structure. Section B 1S a summary of commands available. 
Section C provides the job control language for executing 
the program in both the batch and interactive modes at NPS. 
Section D contains detailed specifications for each avail- 
able command. Pinally, section E gives direction for 
solving larger problems with CAL-NPS. The majority of this 
appendix was originally published as Reference 1. The 
author wishes to express’ appreciation to Professor Wilson 
for permission to use this material. 


A. FORM AND RESTRICTION OF THE LANGUAGE 


CAL-NPS 1s an interpretive language which is designed to 
Manipulate arrays and matrices, to perform standard struc- 
tural analysis operations and to perform conduction heat 
transfer analysis operations. A CAL-NPS program run 
involves the reading of the input deck once and executing 
the commands designated by the operation cards as they arse 
encountered. Looping operations allow a sequence of 
commands to be executed more than once. 

The input deck is composed of operation cards and data 
cards. The data cards directly follow each operation card 
which requires data (see LOOP operation for exception to 
emis) . The operation card contains the name of the opera- 


tion to be executed, names of arrays associated with the 
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Operation and integer constants. Examples of the general 
form of this card are: 

OP,M1,M2,M3,M4,M5,N1,N2,N3,N& comments 

OP, M1,N1,N2 

OP, N1 

OP 
in which OP is the name of the operation to be executed, Mi 
is the name of an array and Ni is an integer. The names of 
OP or Mi are one to eight alphabetic oor numeric characters 


to be selected by the user. The first character of a name 
must be alphabetic. The sequence of terms OP, Mi and Ni 
must be separated by commas. Characters following a blank 


will be printed as comments in the output from the program 
run. 

If an operation attempts to load or generate an array 
which previously existed, the program will delete the array 
before the execution of the operation. A new array need not 


be the same size of the old array which had the same name. 
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Be SUMMARY OF COMMANDS 


ie 


General Commands 





*% indicates a significant change or addtion in CAL-NPS 


START 
S20 P 
NO 
YES 
LABEL 
REA D 
WRITE 
TIME 
SAVE 
RES UME 
LIST 
HEL P 


LOA D 
LOA DI 
ZERO 
PRINT 
DUP 
ADD 
SUB 
MULT 
TRAN 
SCALE 
SOLVE 
DUP SM 
STO SM 
DUP DG 
STO DG 


Initialize for the next problen 


Normal termination 


Temporary suppression of output 


Restores output 


Print comments 


Change logical device for input 


Change 


logical device for output 


Suppress time printout 


Interrupt a problen 


Continue an interrupted problem 


List arrays 


Access 


eneral 


and storage used 
HELP files 


atrix Commands 





Load user defined real matrix 


Load user defined integer matrix 


Create 
Matrix 
Matrix 
Matrix 
Matrix 
Matrix 
Matrix 


null or unitegatrix 
print operation 
duplication 
addition 
subtraction 
multiplication 


transpose 


Multiply a matrix by a scalar 


Solution of linear equations 


Form sub-matrix from large matrix 


Store sub-matrix in large matrix 


Form row matrix from diagonal 


Store row on diagonal 
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MAX - Evaluate row maximums 

NOR M - Evaluate matrix norms 

INV EL - Invert each term in matrix 

SOREL -~ Square root of each term in matrix 

LOG - Natural log of each term in matrix 

PROD - Evaluate product of all terms in a matrix 
DELETE - Delete matrix from storage 


3. Static 


[~~ 


nalysis Operations 


NODES - Input structural joint geometry 
BOUND - Specify boundary conditions 
BEAM - Form 3-D beam stiffness matrix 


TRUSS - Form 3-D truss stiffness matrix 

PLANE ~ Form 3 to 8 node plane stiffness matrix 

SLOPE -~ Form stiffness matrix from slope/deflection eq. 
FRAME -~ Form 2-D frame stiffness matrix 


LOADS - Form load vector 
ADDSF - Form global stiffness and mass matrices 
ADDK ~ Add element matrix to global matrix 


MEMFRC <=- Calculate element forces from joint displacements 
DIS PL 
FORCE 


Print joint displacements 


Evaluate and print member forces 


4. Dynamic Analysis Operations 


FUNG Generate equal interval time function 

STEP ~ Integrate dynamic equilibrium equations 

EIGEN ~ Evaluate mode shapes and frequencies 

DYNAM ~ Evaluate uncoupled equations of motion by 
mode superposition method 


PLOT =“"hane pranter plot of joint time history 


HTXFR - Initiate heat transfer problem 
Coo RD ~ Input nodal coordinates 


ELCON ~ Input element connectivity matrix 
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PROP 
CTE MP 
PROF 
Sie 
USYMC 
LCA P 
CCAP 
FOR M 
CALC 
ODE 
PTEMP 
TOL 
CON V 
DTIM 
ADTIM 
PROMPT 


6. 


GRA PH 
TITLE 
PLH X 
PLST 


LOO P 
NEXT 
Ste. P 


USERA 
USERB 


- Input material property data 

- Input constant temperature node data 

- Establish profile of equations 

- Create symmetric conductance matrix 

- Create unsymmetric conductance matrix 

- Create lumped capacitance matrix 

- Create consistent capacitance matrix 

- Create flux vector 

- Solve time independent systems of equations 
- Solve time dependent systems of equations 
- Print nodal temperatures 

- Set solution convergence tolerance 

- Perform temperature convergence test 


- Set time step increment 


Advance time by one time step 
- Suppress/restore prompts 


erations 


: 


Graphics 


Initiate graphics 
Label mesh plot 
- Plot heat transfer mesh 


- Plot structural mesh 


Loop Operations 





- Start of loop 
- End of loop 
- Conditional skip of operations within loop 


Names Available for User Subroutines 
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C. EXECUTION OF CAL-NPS 


For the time sharing (CMS) system at Naval Postgraduate 
School, do the following: 
_(Use the standard LOSON procedue) 
_link 0040P 191 199 
ENTER PASSWORD: 
_XXXXX (ESAN) 
Rs; T=0.01/0.01 11:09:55 
_access 199 C (Note: You must access the "Cc" 
C (199) R/O disk) 
R; T=0.01/0.01 11310:04 
mead 
ENTER TERMINAL CODE: 


PLOT-10 Compatible Terminal (GRAPHICS) 
IBM 3277 DUAL SCREEN (GRAPHICS) 
Any Alpha Numeric Terminal (NO GRAPHICS) 


_1 (or 2 or 3, as appropriate) 

(The computer responds with several lines of procedure) 
EXECUTION BEGINS ..caece 

ee a SS SS Q.0 SECONDS 
(You are now under the control of CAL-NPS) 

start 

aeaSTART 

(Your own CAL-NPS program is inserted here) 

stop 

"STOP 

R; T=0.01/0.01 11212345 

_log (Terminates session) 
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D. CAL-NPS COMMAND SPECIFICATIONS 


CAL-NPS has most of the standard matrix operations 
plus some special array operations which are useful in engi- 
neering analysis. The following is a list of approximately 
32 operations which are used for control and general matrix 
manipulation. 

A "+" indicates the formation of a new matrix. A 
matrix previously defined with the same name will be 
deleted. A “= yundicates modification of an existing 
matrix. 

Note: Whenever the expression "card" is used it is 
meant to also stand for "instruction" in interactive mode. 
Bee 


as gucrte2ou eliminates all arrays which were previously 
oaded or generated. 


This operation causeS normal termination of a CAL=-NPS 
program. 


NO 
YES 
These orferations are used to selectively suppress output 
from CAL=-NPS. The NO operation suppresses all printing, 


except diagnostics, until the operation YES is encountered. 
Therefore, in subsequent runs o the same CAL-NPS program 
output which was previously correct need not be reprinted ic 
these cards are inserted in the data deck. 


LABEL,N1 


This operation will read_and print N1 comment cards. which 
follow the operation card. Olumn 1 of each card will be 
interpreted as a standard carriage control symbol (i.e. 0 
for double space and 1 f9r skip to the top of the next 


page) . 
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READ, N1 
oe OPERATION IS VALID ONLY WITH THE CP/CMS TIME SHARING 
This operation permits the selection of a user's file or the 


terminal as the input file device. The default is _the 
terminal. If N1 is &, subsequent commands will be read from 
FILE. FTOYFOO01 on the user's A-disk. TEN 2s 5, the 
terminal will be restored as the input file device. All 


disk files prepared for use with this command should end 
with either STOP or READ,S. This command will not be 
executed on the MVS (batch) systen. 


WRITE,N1 


ae 2 OPERATION IS VALID ONLY WITH THE CP/CMS TIME SHARING 


This operation permits the selection of a file on the user's 
A-disk or the terminal as the output file device... The 
default is the terminal and all 2rror messages will (be 
printed at the terminal regardless of the output file device 
selected. If N1 is 8, subsequent pope will be written into 
FILE FTO8FO01 on the user's A-disk. me Nip isi 6, the 
terminal will be restored as the euepas file device. This 
command will not be executed on the MVS (batch) systen. 


TINE 
This operation permits the time printout to be. suppressed 
without loss of other output. A second .TIME will restore 


the time printout unless the print output is suppressed with 
the NO command. 


SAVE or SAVE,M1 or SAVE,M1, N1 


This operation saves all arrays in. storage at. the time of 
1ssuance. Saved arrays will contain all modifications made 
Since their creation. M1 and N1_ are optional and if not 
included, the ee will be stored in FILE 02 on the user's 
A-disk. The saved files will be assumed to be for general 
matrix manipulation or a structural problem. M1 is the file 
name (up to six letters under which the user wishes to 
store the arrays. The file type will be SAVE. Trond 1541, 
a general matrix manipulation or structural problem is being 
saved. If N1 is 2, a heat transfer problem is being saved. 


BeSsUNE or RESUMNE,M1 or RESUME,N1,N1 


This ep ea aon reads a saved file into memory... Any arrays 
currently in storage will be destroyed. A file must have 
been previously created on the user's A-disk using the SAVE 
Operation. M1. and N1 are optional and if not included, FILE 
O02 will be read. It will be assumed to be. a general matrix 
Manipulation or structural problem. M1 ais the file name 
assigned to the saved file on the user's A-disk. The file 
type must be SAVE. If N1 is 1,. a general matrix manipula- 
tion or structural problem is being resumed. PEON A 2S: 2e2a8 


heat transfer problem is being resumed. 


eS T 
The LIST operation prints the directorv information for 
arrays in storage and the amount of storage used. 
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+ 
LOAD, M1,N1,N2,N3 


This operation will jload an_array of real numbers named M1 
which has N1 rows and N2 columns. The terms of the array 
are punched, in row-wise sequence on data cards following 
this operation. N3 is optional. 1ivweNSs 25 Zero or Diank, 
the cards_are punched in(8P10.0) format. If N3 is one, an 
additional card which contains the format of the data cards 
must preceed the data, For example, if the data is to be 4 
numbers per card in field widths of 15, 
would contain the ee ne information: (4F15.0). If N3 is 
nine, the data cards will be read in free format. 


+ + 
LOADI ,M1,N1,N2,N3,N4 or LOADI,M1,N1,N2,,N4 or 
LOADI,M1,N1,N2,N3 or LOADI,M1,N1,N2 or 
LOADI,M1,42,N1,N2,N3, N4 : 


This operation will load an. integer array named M1 which has 
Nt crows and N2 columns. |The tefms of the array are punched 
in row wise sequence on data cards following this operation. 
M2,N3 and N4Y are optional. If N3 is zero. or blank, the data 
must be punched in, (1615) format. If N3 is one, an addi- 
tional card containing the format of the data cards must 
follow this operation and preceed the data. For example, if 
the data is to be 4 numbers er card in field widths of 10, 
the additional card would contain: (4110). If N3 is nine, 
the data will be read in free format. LD ss ; 
me the letters "ENCR™ are placed in the position of M2, this 
operation has an increment generation capability. Data must 
be entered as follows: 
Item Contents 

1 Row number 

2 Value 1 

3 Value 2 


etc. 

N2+1 Value N2 

N2+2 Generation code 
If the eneration code is not zero, the next card must 
contain the following: 

Item Contents 
Row number increment 
2 Value 1 increment 
3 Value 2 increment 


etc. 
N2+1 Value N2 increment 
N2+2 Last row to be generated 


+ 
ZERO, 41,N1,N2,N3,N4 


A real matrix named, M1 is created with N1 rows and .N2 
columns. The terms in this matrix will have the following 
values: 

aa eee ay = N3 Ee= 16 se 

M1(I,J) = N& J = 1,... 
Therefore this operation can be usei to form null or unit 
matrices. 
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PRINT,M1 or PRINT, M1, N1 or PRINT,M1,N1,N2 or 
PRINT,M1,N1,N2,N3 


This operation will print the array named M1 in a matrix 
format of as to elght columns per _ line. N1,N2,N3 and N& 
are optional. N1-is the number of comment cards (followin 
the operation card) which will _be read and printed. N 
defaults to zero.. If N2 is included, the matrix will be 
pee ce, 1a partitioned form with N2 columns per partition. 
ines will -have N2%15 + 5 characters. N2 defaults to 8, 
printing 125 charcters per line. If N3 is included and is 
greater than zero, integer format ee is used. The default 
value is zero and real format (PD15.7) 1s used. The user is 
cautioned not to overcome the eobact ey of the displaying 
device in use to avoid wrap around on the screen. 


+ 
DUP ,M1,M2 


This operation will form an array named M2 which is iden- 
tical to the array named M1. 


ADD ,M1,M2 


This operation will replace matrix M1 with the sum of the 
matrices M1 and M2. 


SUB,M41,M2 


This operation will replace matrix M1 with matrix M1 1éess 
Matrix M2. 


+ 
MOLT, M1,M2,M3 


This operation, generate 
an 


x M3 which is the 
product of matrices M M2 


= ae ap 2 2 2 ow eee 2 Ce «ee a 2 eo 6 eee ee ee ew OP ee ee 2 ee ee ee ae 6 ee 6 ee a ew DD ow aD om 2 em oD aD ew 6 ow a ep em ep 6 eo op eo a a am @® a a= 2D 


+ 
TRAN, M1, M2 


This operation generates a new matrix M2 which is the trans- 
pose o£ matrix M1. 


SCALE,M1,M2 


This operation replaces each term in the matrix named M1 
ee, RS term multiplied by the term M2(1,1) of the matrix 
name s 


42 





SOLVE,M1,M2,N1,N2 or SOLVE, M1,M2,,N2 or SOLVE,M1,N1,N2 or 

SOLVE ,M1,42,N1 or SOLVE, M1,M2 

Tf N1 = 0, this operation solves the matrix equation AX=B. 

M1 1s the name of the A matrix and M2 1s the name of the B 

Matrix. Matrix A is triangularized and the results, X, are 

stored in M2. If N1i=1, Matrix A 1s triangularized only. 

Tf Ni=2 for a oe B matrix and the Matrix previously 

trianguiarized, the B matrix is replaced by the results, X. 
N atrix A is replaced by its inverse FOR SYMMETRIC 


=3 M 

MATRICES ONLY. 
If N2=0 or blank, matrix A is symmetric. If N2 is nonzero 
the matrix A is not symmetric. ; 
For symmetric matrices, A is. factored into the LDL forn. 
The diagonal D matrix is stored on the diagonal of A. The 
parameter N2 permits the direct solution. of non-sSymmetric 
Systems of equations. . If N2 is not equal to zero, an LU 
decomposition of matrix A. will be performed. No direct 
replacement of M1 by its inverse is available for the non- 
Symmetric case. . Instead use the ZERO operation to create 
an goguT tty Matrix M2 of the same order as M1. The command 
SOLVE,M1,M N2 will then replace the matrix M2 with the 
inverse of the matrix A. 


+ 
DUPSM,M1,M2,Ni1,N2,N3, NU 


This operation forms a new submatrix named M2 ,with N3_rows 
and NY columns from the terms within the matrix named M1. 
The first term of matrix M2, M2(1 ty will be from row N1 
and column N2 of matrix M1, or M1 Ni, 2). 


Smests41,M2,N1,N2 
This operation stores a submatrix named M2 within the matrix 
named 41. The first term_of the grup aerik M2 will be stored 


at row N1 and column N2 of matrix M1. The terms within the 
area of M1 in which M2 1s stored will be destroyed. 


+ 
DUPDG,M1,M2 


This operation forms a new matrix named M2 from the diagonal 
terms of M1. 


STODG ,M1,M2 


This cperation, stores a row or column matrix named M2 at the 
diagonal locations of matrix M1. 


This operation forms .a column matrix. named M2 in which each 
row contains the maximum absolute value of the corresponding 
fOw in matrix M1. The maximum and its column number is 
printed for each row. 
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+ 
NORM, M1,M2,N1 


If N1 = O, .a row matrix. named M2 is formed in which each 
column contains the sum of the absolute values of the corre- 
sponding column of matrix M1. If NI is not equal to 0, a 
row matrix named M2 is formed in which each column contains 
the square root of the sum of the squares of the values of 
the corresponding columns of matrix 41. 


SQREL,M1 


This operation replaces each term in the matrix named M1 
with the square root of the term. 


This operation replaces each term in the matrix named M1 
ee the natural tba Of the term. 


_ Eee ape eae ee ene ee ae eee SP ae ee ee Pee SP eee Se SS ee PSP SPS SP SPP PS SP SBS SSP SESE SSP SPS SPS SSP SS eS Se ee ee ee es 2 oe Se aD 


+ 
PROD, M1,M2 


This operation forms a1x 2 array named M2 which contains 
the product. of all.terms in the matrix named M1. The 
product, X, is stored as two numbers of the form: 

= p10arrR 
in which M2(1) = P and M2(2) = E, the exponent. 


DELETE, M1 


This operation. will cause the elimination from storage of 
the array named M1. 


CMS ,N1 

This operation allows the user to issue CMS commands while 
under the control of CAL-NPS. In general a command that 
reloads the virtual core will not be allowed. Examples are 


PORTHX, SCRIPT, XEDIT, any pore Ae Cee SOL: SORT, and 
other farge system modules such as COPYFILE and MOVFILE. N1 
is the number of words in the command, (1 *o 9 are allowed). 
Note that a parenthesis - "(" or ")" is counted as a word. 
All words must be left justified. 


The user will be pron ted for each word. If this operation 
ae used in an FTOQFO0O1 FILE, each word must be on a Separate 
ine. 
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2.  ot@etlcuknatysis Operatzons 


The purpose of this series of operations is to form 
the total stiffness and diagonal (lumped) mass matrices for 
systems of two- or three-dimensional elements. For three- 
dimensional analysis there are beam and truss elements 
available. For two-dimensional analysis, there is a frame 
element, a slope/deflection element for beams, anda &4& to 8 
node isoparametric finite element available. 

After the creation of an array containing the coor- 
dinates of the joints of the systen, the specification of 
displacement boundary conditions, the tabulation of material 
and section properties, the mass and stiffness matrices are 
formed for each structural member and placed in sequence on 
low speed storage along with the global equation numbers 
which are associated with their stiffness terms. In addi- 
tion, the member force-dis placement transformation matrices 
are formed and stored on a separate low speed storage file 
along with the appropriate displacement numbers. 

The NODES operation is used to specify or generate 
the geometry of the systen. The operation BOUND specifies 
which joint displacements exist and assigns internal equa- 

ion numbers to these displacements. Therefore, each joint 
may have from zero to six displacement degrees of freedon. 
Tables of material and section properties for the various 
members are loaded and printed as standard arrays of 
information. 

A special operation, ADDSF, is used for the direct 
addition of element stiffnesses to form the total stiffness 
and diagonal mass matrix of the systen. The ADDK operation 
may be used to add individual elements into the total systen 
matrices. The LOADS operation specifies the concentrated 
joint loads for all load conditions. After the direct solu- 


tion for joint displacements due to static or dynamic loads, 
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the member forces can be evaluated using the MEMFRC opera- 


ON . The DISPL operation is used to print the displace- 
ments in joint number order. 
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* “ 
NODES,M1,N1 or NODES, M1,N1, N2 


This operation generates an array (N1,3) named M1 which 
contains the coordinates for all joints ina structural 
system. N2 1S optional. Data must be entered in free 
format as follows: 
Item Contents 
1 Node number 
2 X-coordinate 
3 Y-coordinate 
4 Z-coordainea <= 
If N2 = 1, there is a joint generation and coordinate systen 
conversion capability. Data mu be entered in free format 
as follows: 
Item Contents 
1 Node number 
2 X-coordinate 
3 Y-coordinate 
uy Z-coordinate 
5 System type 
6 Generation code 
System ae refers to the oo used when inputing the 
data. A coordinates will be converted to ‘the cartesian 
system for use by CAL-NPS. 
system type Systen, 
1 Cartesian ; 
Zz Cylindrical, Z axis longitudinal 
3 Cylemere cal, toax.s Longitudinal 
4 Cylindrical, X axis longitudinal 
5 Spherical 
The input data is the same as above with the following 
correspondence: 
Cartesian Cylindrical Spherical 
x r 0 
Y Q 0 
Z Z d 
Tf the generation code is not zero, the next card is.a 
generation vector for the self generation of nodes. = 2s 


ormatted as follows: 


item Contents ; 
1 Node number increment 
Z X increment 
3 Y increment 
4 2 increment 
5 Last node number to be generated. 


It is assumed that the increments 
of reference as the preceeding car 
This operation must be terminated by a 
zeros and blanks. 


aoeooo to the same system 
line of alternating 
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| + 
BOUND,M1 
This operation specifies the displacements which are nonzero 
for the structural system of joints specified by the NODES 
operation. Where: 


M1 = Name of boundary condition code array to be 
generated. 


This ee ee followed ey a series of cards containing 
the following information in free format: 


Iten Contents 
1 Node number for the first node in a 
series of nodes with ildentica 
displacement specification. 


Z Node number fcr the last node in the 
series. 


X-translation 
Y-translation 
Z~-translation 
X-rotation 
Y~rotation 
Z-rotation 


Oo ON KH MN F&F W 


Node. number increment u 


sed ¢ 
conditions for addtional n 


-O generate 
odes. 
ero or unde- 


A translation or rotation equals: (a) zero for z 
displacements to 


fined pees ee or (b) one for nonzero 
be evaluated by other operations. 


If a node boundary condition is not speeltved. all displace- 


ments at that node are assumed zero. Cards may be supplied 
in 28 order. If node boundary conditions are specifie 
more a 


nh once, the last definition is used. This sequence 
of data must be terminated by aocard of alternating zeros 
and blanks. 


The selection by the aa of which nodes have nonzero 
isplacments requires an understanding of the direct stirf- 
ness precedure. Displacement degrees of freedom which have 
no stiffness associated with ‘the displacement must be 
considered to be undefined since it iS not possible to 
develop an equilibrium equation for that direction. The 
total number of nonzero displacements specified will be the 
size of the total stiffness matrix to be defined by the 
ADDSF operation. 
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+ 
BEAM, M1,N2,M3,M4 


This operation calculates the element stiffness, mass and 
force-displacement transformation matrices for 3-D beam 
members. These arrays are stored in sequence on low speed 


storage to be used by other operations where: 


is ¢t name o e coordinate .art 

M3 is the name of the boundary condition array 

M4 is the name of the array which contains beam 
properties and has, been loaded by the standard 
matrix LOAD operation 


M1 +8 pe nane of Epe beam element oe e 


oy card for each, beam in this grou Qf beam elements must 
trollow this operation. The bea cards are punched in free 
format, wheres: 


Item Contents ._._. ; 

Beam identification number 
Node number I 

Node number J 

Node number K 

Beam property number NP 


‘~< 


nN uN) 


This sequence of cands must be terminated with a card of 
alternating zeros and blanks. 


The material and geometric properties for each element are 
given in the MU array in the following order: 


M4 (NP,1) = Axial area of member, A . 

MG (NP,2) = Torsional moment of intertia, J 
M4 (NP, 3) = Moment of inertia about axis 2, I 
M4 (NP,4) = Moment of Lnertia about axis 3, I 
M4 (NP,5) = Modulus of elasticity, £ 

M4 Ree = Shear modulus, 


G 
Mass per unit length of bean 
where NP is the specific material property number specified 
2 


in item 5 of the beam card. an ocal Sign convention is 
given in the following figure. 
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LOCAL BEAM REFERENCE SYSTEM: 
AXES | AND 2 ARE IN THE 
PLANE DEFINED BY PLANE 
[,J ANO K NODES. AXIS 1 
[S DEFINED BY LINE I-J. 
AXIS 3 IS PERPENDICULAR 
TO THE [, J, K PLANE. 


ry 
: 7 K 
é \ i 
B\ RL ee 

F Pe 

4 as 

= 

\ F 
Fg 


DEFINITION OF POSITIVE BEAM FORCES 
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+ 
TRUSS 81 ,N2,43, M4 


Tbis operation forms the element, stiffness mass and force- 
isplacemen* transformation matrices for 3-D truss members. 
The arrays are stored on low speed storage in sequence and 
will be used by other structural operations. 


1 is «he name of this group of truss members 

2 iS the name of the céordinatse array 

3 is the name of the boundary condition array 

G& is an NP by 3 array of section Be ears a 
which NP is the number of different section 
properties and 


Se 
ne 


M4 (NP, 1) = The cross-sectional area, A 

2) = The modulus of elasticity, £E 
_ M4(NP,3) = the mass per unit length of *he member. 
This matrix can be loaded by the matrix LOAD operation. 


his operation is followed b ne card,per truss member in 
ree format with the fal vewine Vnformati tie 
item Contents ae ; 
Truss member identification number 
2 Joint number I 
3 Joint number J 
4 Section property number, NP 


This operation must be terminated by a card of alternating 
zeros and blanks. 


+ + 
FRAMNE,M1,M2 
This eration orms the 6x 6 stiffness matrix for the 
aoa hennionel foume member shown Bade’ 
Y 


F 


YJ 





FRAME MEMBER 
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The properties of the member are defined on one cards imme- 
diately following, the FRAME operation card. This second 
card is punched in free format and contains the following 
information: . 


Items Cone nts 






1 Xial area, A ~ a 
Z Modulus of elasticity, E 
3 Moment of inertia, I 
a Xi 
5 Y 
3 i} 
7 Y 
Y 
Eee 
es 
Uys 
Pcnmeaniecwec ye 
ra e 





x 


Xy xs 


GEOMETRY AND JOINT OISPLACEMENTS 


M2 is a 3x 6 force-displacement transformation matrix which 
a sosec. on the positive definition of the element forces 
shown below. 


M P 


J 
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These forces can be calculated from the following matrix 
equation with the MEMFRC operation. 


Uxi 
Uyi 
= [m2]] Us 
[ Mi 
Uy 
Uei 


ee are, ab SDP SD aD 6D SPD ODP 24D 22a ae 4 aE ae SGD apa a GD SP GD SPE SEB BE EE Bee ee eee ae SE RPE ew 8 BeBe Pe OS ae SSeS eS 


wis 


+ 
SLOPE ,M1 


This operation cores aux & stiffness gatrixe M1 for a bean 
or column member from the classical slope-deflection equa- 
tions. The poe 9 of the member are defined on one Card 
immediately Following the cperation. This second card is 
unched in free format and contains the following informa- 
ion: 
Item Contents. ; 
Moment of inertia, I 
ze Modulus of elasticity, £ 
3 Length of member, L 


The sign convention is defined as follows: 





[I 











‘T 
:% 
T i 
DISPLACEMENTS 
f s 
‘ ———— V2 
FORCES FORCES OISPLACEMENTS 


The member forces are defined in terms of joint displace- 
ments by the following slope deflection equations. 


i Et |40,+20-$( ur) | 


Mz = Ellog,+40,- 74-2) 


Vy = ~V_= Ms+M2 
L 
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* 
PLANE,M1,M2,M3, M4,N1, N2 


This operation calculates the, element stiffness,mass and 
stress-displacement transformation matrices for 4 to 8 node 
isoparametric elements (Y-Z plane oe -These arrays are 
stored on low speed storage *o be used later by other, opera- 
tions (1.e., ADDSF and FORCE). The arguments are defined as 


M1 is the user defined name of the slement group 

M2 1s the name of the joint coordinate array 

M3 is the name of the Boundary condition array 

M4 is the name of the array which contains the 
material properties of the elements (one row per 
different material) where 


M4 (NP, 1 
M4 (NP,2 
M4 (NP, 3 


M4 (NP,4 


Modulus of elasticity, E 
Poissons ratio, Vv 
Thickness of element 


Mass density of the element. 


NP is the material identification number. 


Ni and N2 are the number of integration points in the r and 
s directions respectively. 


eee card for each 3 to 8 node element in the BINS must 
follow the operation card. The cards are punchéd 1h free 
format and contain the follecwing information: 


Item Contents. ae ; 
Element identification number 
Node number N1 
Node number N2 
Node number N3 
Node number N4 
Node number N5 
Node number N6 
Node number N7 
Node number NS. 
Material identification number, NP 

Natural Coordinate of stress output r1 
Natural Coordinate of stress output sl 
Natural Coordinate of stress output r2 
Natural Coordinate of stress output s2 
Natural Coordinate of stress output r3 
Natural Coordinate of stress output s3 


ANSE WN 32 OW O~YHNUNE WN 


N4 through N8 are optional, but zeros must be inserted for 
them if unused. he midside nodes, if present, must be 
within the center half of the side. The local numbering 
System for the element is shown in the following figure. 
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ISOPARAMETRIC ELEMENT 


Stresses will be printed aye ere FORCE operation at the three 


Olnts defined in items 


through 16. The forces are 
efined as follows: 


1 
F, Ox x 
1 

F. Oyy 
F, T xy 
2 

F, - OF» x 
= 2 

Fs Oyy 
Fe Txy 
3 

FY Ox X 
3 

Fe Ovy 
F, Try 
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‘ 
LOADS ,M1,42,N1 


This operaticn ,forms a load matrix named M1 of N1,. columns 
(N17 load conditions) where M2 is the name of the boundary 
condition array generated by the operation BOUND. This 
operation is followed 23 a2 series of cards - one for each 
loaded joint for each load condition. These cards are 
punched in free format as follows: 


Item Column 
Joint number 

Load conditior number 
Load in X-direction 
Load in Y-dir ection 
Load in Z-direction 
Momemt about X-axis 
Momemt about Y-axis 
Momemt about Z-axis 


OYAnANIE Wh 


This series of cards must be terminated by row of alter- 
hating zeroes and blanks. 


+ + + 
ADDSF,M1 or ADDSF,M1, M2 


This operation forms the total stiffness matrix named M1 and 
the lumped mass matrix name for the structural systen 
from the element stiffness and mass matrices which ar? 
stored on low speed storage. These matrices can be printed 
with the PRINT operation. Peeiceis NOt Spectficd, he row 
mass matrix M2 will net be formed. 


= eeu ae ee eae ee aan 2 eee eee awa ae 2a See Pee PSPS PPP SP PSs SP PPP PP SP SPS Pee SS SP Oe OS See ae ee a a = ae = a= 


ADDK,M1,M2,M3,N1 


This operation adds the element stiffness matrix named M2 to 
the total stiffness matrix named M1, wher M1 was previously 
defined and initially set to zero. M3 is the name of the 
integer array in which the column number Ni contains the row 
or column numbers in the total stiffness matrix where the 
element stiffness terms are to be added. 


DISPL,M1,M2 


This operation prints the displacement vector named M1 in 
jOint sequence order, where M2 1s the name of the boundary 
condition array. 
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+ 
MEMFRC,M1,M2,M3,M4,N1 


This operation eee the element stiffness matrix, named 
M1 by ‘the joint displacement matrix named M2. M3". 21S the 
Mame of the integer array in which the column number N1 
contains the row numbers,in the displacement matrix, M2, 
which are to be multiplied by the element stiffness (or 
force-displacement) matrix, M1. The results of this nulti- 
plication are stored in the array named M4. 


+ 
FORCE ,M1,M2,M3 or FORCE, M1, M2 


This operation,calculates the member forces for a group of 
elements in which 


M1 is the name of she element group 

M2 is the name of the displacefent matrix 

M3 is the name of the matrix_in which the forces are 
stored in the order calculated. 


If M3 is not specified, the element forces will be print 
only and will not be retained in storage. For t i 
element only the member axial force, fF, will be calculated 
for each member. For the BEAM element, eight forces will be 
rinted with reference to the positive definition shown in 
he BEAM operation. 
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3. Dynami 


nalysis Operations 


(i 


_The following operations were designed to evaluate 
the dynamic response of structures subjected to arbitrary 
time-dependent loads. If these operations are used in 
connection with the standard matrix operations and the 
structural analysis operations, a dynamic analysis is a 
relatively simple procedure. The user has the option of 
using the mode superposition method or a direct step-by-step 
integration of the dynamic equations of metion. The user 
may examine the spectra of both input loading and calculated 
displacements. In addition, the contributions of the indi- 
vidual modes may be evaluated and compared. 

The most common and convenient form for time- 
dependent data to be specified is as straight line segments 
between given time points. Therefore, an operation which 
generates values at equal intervals is necessary. Another 
common characteristic of time-varying loads on structures is 
that it is normally possible to represent the loads at all 
points on the structure by the product of two matrices, a 
column matrix indicating the spatial distribution of loads 
times a row matrix which indicates the values as a function 
of various times. If a more complicated loading is 
required, it is possible to perform more analyses, each 
Within the restrictions of the program, then add the results 
of each analysis. 

The following operations have been added for the 
major purpose of perfcrming dynamic analysis. 
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+ 
FUNG, 41,82,43,N1, N2 
This operation nerates 2 matrix named M2 which contains 
values at egual TPatervais of the function apendied in the 


peel named M1. The array M1 must be a 2 by K array of the 
orm: 


M1] = 
ie er om ff cee £, 


ach numerically represents a function of the form shown 
elow: 


Pel ) 
Fal As 


rar 
FAT, 
wie : 


V V > | 


The time interval t is specified in the 1 by 1 matrix named 
M3. N1 specifies the total number of values. to be gener- 
ated, and is_the number of columns in M2, If N2 = 0, _.the 
array 42 will be atx, N11 row matrix =in which the first 
value will be ff. If N2 is not equal to 0, the array M2 will 
be a 2 x N1 matrix of the following forn: 


€, t, + At t, + 2At eoecece 
f,  £(t, + At) £(t, + 2At) ..... 


[m2] = 
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- - + 

STEP, M1,M2,M3,M4,N5,M6,M7,4N8,N1,N2 

This operation calculates the dynamic respcense of a 
s 


St 
tural system using direct step-by-step integration of 
following linear matrix equation of motion: 


CM){U} + (CJ(U} + CKJCU} = R(t) = PF(t) 


ruc- 
the 


Wheres, . 
is the name of the N x N stiffness matrix K 
M2 is the name of the N x N mass matrix 4M 
M3 is the name of the N x N damping matrix Cc : 
M4 is the name of the N x 3 initial condition matrix 


0 in which: . 
I,1) 1s a vector of displacements JU 
cae: is a vector of velocities JU 
(1,3) is a vector of accelerations JU 


MS is the name of,the N by N2 matrix of calculated 
displacements in which column 1 represents the 
displacements at time i1°N1* t | 

M6 is the name of the N x 1 load distribution matrix 


Aaa 


P 

M7 is the name of the 1 x k row matrix representing 
the load Boe eS at equal time increments 
?, where k = N2/N1 . . 

M8 is the name of the 1x 1 matrix containing t 

N1 is the output interval for the displacements 

N2 is the total number of displacement vectors to 
be calculated. 


The total time for which results will be calculated by this 
Operation 1S NI¥*¥N2" t+. This operation must be followed with 
Bee cate card in free format containing the following infor- 
mation: 


Iten Contents 
1 DELTA 
2 ALPHA 
3 THETA 
Different values of delta, alpha and theta wil 


1 allow the 
user to select different methods of BES ovac coe saa 
tion. The following table lists some possibilities: 


DELTA ALPHA THETA 
Newmarks Average Acceleraticn 1/2 1/4 1.0 
Linear Acceleration 1/2 1/6 V0 
Wilson's Theta Method ee damp ing) 12 1/6 i= 
Wilson's Theta Method (high damping) we 1/6 ZO 
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—_ + ae 
EIGEN,M1,42,M3,N1 


This operation solves the following eigenvalue problen: 
K 6 = M®@X 


on which the N x N, symetric, posj;ive-semideft-nits Matrix K 
is named M1. fhe matrix M is a fagona matrix of nonzero, 
positive terms designated by M3. he matrix M3 must bea 
row or column matrix containing only the diagonal terms of 
Me. The eigenvalues, , are stored in matrix M3. The eigen- 
values are ordered in Dene eee SY increasing order and the 
elgenvectors » are stored in the corresponding columns of 
the matrix M2. _ The number N1 specifies the “approximate 
number of Se ene figures of the eigenvalues. Lf Ni 2s 
zero or blank, 4 figure accuracy will be used. The maximum 
accuracy possible 1S 16 figures. The use of more than 12 
figure accuracy is not recommended. 


The Bode es reduces the problem to standard eigenvalue forn 
by the following transforma tion 


where 
Ke =n'K o 
I = mn’ M m 
in which 
M, = 1/ v Mii 


The calculated mode shapes, , are normalized as follows: 
@MO=T o'™™* 6 = 


The program uses the standard Jacobi diagonalization method 
to solve for all eigenvalues and sigenvectors. 
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+ 
DYNAM ,M1,M2,M3,M4,M5, M6, N1 


his .operati evaluate the following. set of uncoupled 
second rder Bcferenttat equations associated with the hode 
superposition method for the dynamic analysis of a struc- 
tural systen. 


%,+ 2A X% + 4X, = RG) i = 1 to N nodes” 


M1 is the name of a row or column matrix which contains the 
N terms (frequencies in rad/sec). M2 is the name of a row 
or column matrix which contains the N terms (ratio of 
modal damping to critical damping). 


‘ie generalized time-varying forces P (t) are not_ specified 
ilrectly put are evaluated” from more fundamental informa- 
UO Nie The forces for all modes. are evaluated at specific 
times by the program from the following matrix equation: 


p= pf = M374 


In which p is a specified N_x 1. vector named M3, and fisa 
1x N row matrix which will be generated from the 2 by k 
array named M4, The array M4 is the same form as the input 
array described under the operation FUNG. t 1s not neces- 
sary to use FUNG before the DYNAM operation. 


M5 1s the name of the N x N1 array which contains the gener- 
alized displacement X(t). 


M6 is the name.of the,1. x. 1 array which .contains the time 
increment associated with the generalized displacements. 


N1 is the number of displacements to be generated. 


The methced of integration used is exact for straight line 
segments. 


PLOT, M1,N1 


This operation will prepare a printer plot of selective rows 
of the matrix named M1. NT 1s the number of rows of M1 
which will be plotted by this operation. This operation is 
followed by N1 cards in £— (1a1,I4%) format with the following 
information: 


Columns Contain 
1 Plot symbol - any peeeunen symbol 
a eS Row number to be plotted 


The aes automatically searches the information to be 
plotted for the maximum and minimum values. The difference 
in these. numbers divided by 120 spaces is selected as «he 
plot scale. 
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G4. Heat Transfer Operations 


The purpose of this series of operations is to form 
the total conductivity and heat capacity matrices for 
systems of two- or three-dimensional elements, to form the 
flux vector and solve the defined set of equations. For two 
dimensional elements there are 4 to 9 node isoparametric 
elements. For three dimensional elements there are 8 to 21 
node isoparamtric elements. 

After the creation of an array containing the coor- 
dinates of the system nodes, the specification of element 
connectivity, the specification of material properties and 
the specification of constant temperature nodes, the conduc- 
tivity and capacity matrices and the flux vector are formed. 
The equations are solved by the appropriate equation solver 
and the temperatures are printed in node order. 

The HTXFR operation initializes the problen. The 
operation COORD is used to specify or generate nodal coordi- 
nates. Element connectivity is specified by the ELCON oper- 
ation. Material properties and element boundary conditions 
are input via the PROP operation. The CTEMP operation 
establishes designated nodes as having constant tempera- 
tures. The equation profile for the problem is generated by 
the PROF operation. 

To form the conductivity matrix, the operation SYMC 
is used. The heat Capacity matrix can be approximated with 
either a consistant (CCAP) or lumped (LCAP) matrix formula- 
tion. The flux vector is formed with the FORM operation. 

Time independent systems of equations are solved 
With the CALC operation. Systems of equations involving «the 
first derivative of temperature with respect to time are 
solved with the ODE series of operations. Nodal temepra- 
tures are printed in node number order with the PTEMP opera- 
ON» 
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his operation initializes the heat transfer problem. he 
ollowing iat omnee ion. entered tn free format must omen 
this operation: 


Iten Contents 
1 Number of nodes (NUMNDP 
2 Number of elements (NUML) 
3 Number of material sets (NUMMAT 
rT Spatial dimension (2 or 3) (NDM 
5 Number of unknowns per node (NDF) 
6 Maximum number of nodes per element (NEN) 


The number of unknowns per node will always be 1 for heat 
transfer probiems. 


This operation creates. an array which contains the coordj- 
nates for all nodes ina _heat transfer system. Data is 
entered in free format as follows: 


Item Contents 

1 Node number 
X-coordinate 
Y-coordinate 
Z-coordinate (if 3-D, else omit) 
System type 
Generation code 


ONS Wh 


System type refers to the coordinate system used when input- 
eomg data. All coordinates are converted to the cartesian 
coordinate system for use by CAL-NPS. 


eee TR Systen 


0 Cartesian . . 
> Gylanedcical, Z-akis longitudinal 
3 Gyiemostcal, Y-ax1s longitudinal 
4 Cylindrical, X-axis longitudinal 
5 Spherical 


The input data is the same as shown above with the following 
correspondance: 


oe oe Cylindrical Spherical 

Y 6 

Z Z, d 
If the generation code is not zero_ (0) 
g 


eneration vector for automatic node 
ree format as follows: 


jee ne Next. cava 15.28 
eneration. It is in 


Iten Contents 

Node number increment 

X increnent 

Y increment. . 

Z increment (if 3-D, else omit) 
Last node number to be generated 


MNG&Wwh— 


This operation must be terminated by a card with alternating 
zeroes and blanks. 
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This operation _creates an array which gontains the element 
connectivit data or the ments in a eat transfer 
system. Data is entered in free format as follows: 


Iter Contents 
Element number 
Node 1 number 
Node 2 number 
etc 
Node N number 
Material set number 
N+ 3 Generation code 


SEQ 
+? ftUWNe 
N33 


If the generation code is not zero, (0), the next card is a 
generation vector o the automatic generation of element 
connectivity. It is entered in free format as follows: 


Iten Contents . 
1 Element number increment 


2 Node 1 increment 
3 Node 2 increment 
etc etc 


N+ 1 Node N increment 
N#2 Material set increment (usually 0) 
N+ 3 Last element number to be generated 


This operation must be terminated by a.row of alternating 
zeroes and blanks. The node numbering conventions fof 
element connectivity are shown below. 





; 





2-D Element 3-D Element 
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This operation inputs the material BREE data for. .a heat 
transfer system. The ember .pLOrme lon, entered in free 
format, must follow this operation: 


Item Contents 
1 Material set number 


2 Element type number (2 for 2-D, 3 for 3-D) 
Additional_information must be provided depending on the 
type of element being used. . ote tha Material property 
data must be in consistent units as shown below: 

eee 5S Units SI Units 
k Bie cee F W7-°7C 
C BTU/1bm-OPF kJI/kg-9C 
rho lbm/ft3 oe 
qtv* “SReyhr-cts /n3 


2-D Elements 
This information must follow the PROP operation card and its 


two required entries. An entry must be made for each iten. 
ee Boe item is temperature dependent, the entry will be 
ignored. 


Item Contents . . . ; 
Conductivity in the X-direction 
Conductivity in the Y-direction 
Specific Heat 
Density . 
Heat generation per unit, volume 
Number of eed re EL Cn points per direction 
se to 6, default is 4 
eometry type (see below) oa 
Total number of lines with specified boundar 
conditions in elements with the same materia 
set numer (NLBC) (see below) 
9 Temperature dependence code (see below) 


O~ BANEWN- 


Geometry type is 1 for plane geometry and 2 for axisymmetry. 


The temperature dependence code is 0 if all material proper: 
ties are constant and 11f any property is temperature 
dependent. If the code 1s 1, the following information (in 
free format) is required. 


SC Contents 


Conductivity in the K=-direction code 
2 Gene tet ayy in the Y-direction code 
3 Heat capacity (specific heat*density) code 
4 Heat generated per unit volume code 


where 0 means a constant property and 1_means a temperature 
meee ces t Bese Temperature dependent properties are 
entered in the form of a. table. The tables are consecu- 
avery 2nput for conductivity in the X-direction, conduc- 
feyity in the Y-direction heat capacity and heat generated 
bee uhit volume. Omit the tables or which the temperature 
ependence code is zero. Tables are input in free format as 
shown on the following page. 
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[tea Contents 

Number of data pairs to be entered (This 
should be a single card) 

Temperature 1 as 

Property 1 {These two entries should be on 
one card) : 

Temperature 2 

Property 2 (These two entries should be on 
one card) 


c 
Temperature N ; 
1 Property N (These two entries should be on 


one card) 


n ae Be havea egitied oundary condition (NLBC>0),. 
aa be subaitted ak eet y rg Paco se se Tae de 
car 4 bust to more than ae “boundary, condition, a card must 
be used for each one of these conditions. The total number 
of cards must equal NLBC. The information is entered in 
free format as follows: 

Iten Contents 
Element number 
he gels ES eo ares code (see below) 
Line code (see below 
Property value (see below) 
Ambient temperature 
The beundary ao codes are: 


By 0 
Zeit NE Who = 
$20 


2 


WN WIth) = 


1 ux 
2 Catenion (constant coefficient) 
: Radiation 


Convection (temperature dependent property) 
The line poet ares 


1 =z +1 line 
=-1 = -1 line 
2 TT = #1 Line 
-2 T = -1 line 


Aa eS 


(O-----b 4 





Qa 


The property values are; 
Flux - Flux per unit area 
Convection - constant heat transfer coefficient 
ignored if temperature Sous) 
Radiation - Product of emissivity by Stephan- 
Bolt zman constant 
ct temperature is ignored for the flux boundary 
conga. ti 
If the boundary condition code is 4 (temperature IRI 
heat transfer eo eeece ent) a table must follow for the 
temperature dependence. f is input in free format as is 
shown on the following page. 
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Item Contents ; : 
1 Number cf data pairs to be entered (This 
should be a single card) 
Temperature 1 ele : 
Heat transfer coeficient 1 (These two entries 
should be on one card) 
Temperature 2 , 
Heat transfer coeficient 2 (These two entries 
should be on one card) 
etc 
Temperature N aan ‘ 

1 Heat transfer coeficient N (These two entries 
should be on one card) 


waar WME Wr 


Srow 
+20 


3-D Elements 


This information must follow the PROP op 


eratzon card and its 
two required entries. An entry must be made for each iten. 
If the item is temperature dependent, the entry will be 
ignored. 

Item Contents , . . . 

1 Conductivity in the X-direction 

2 Goneuctivicy in the Y-direction 

3 Conductivity in the Z-direction 

4 Specific Heat 

2 Density . 

6 Heat generation per unit. volume : 

7 Number of integration points per direction 
A to 6, default is 4 

8 eometry type AP Se below) | a. 

9 Total number of surfaces with specified 
eee conditions in elements with the same 
material set numer (NSBC) (see below) 

10 Temperature dependence code (see below) 


Geometry type is 1 for plane geometry and 2 for axisymmetry. 


The temperature dependence code is 0 if all material proper- 
ties are constant and 1 i1£ any .property 1s temperature 
dependent. If the code is 1, the following information (in 
free format) 1S required. 


Item Contents . . , ; 

Comauct tvicy anesthe X=d*section code 
Conductivity in the Y-direction code 
Conductivity in the Z-direction code 

Heat capacity (specific heat"density) code 
Heat generated per unit volume code 


(ni WN 


where 0 means a constant Be eee ane tT means a temperature 
Bee oocent Cee es Temperature dependent properties are 
entered in the form of a. table. The tables are consecu- 
Pevely input for conductivity in the X-direction, conduc- 
meyrty in the Y-direction, conductivity in the Z-direction, 
heat capacity and heat generated per unit vclume. Omit the 
tables for which the temperature dependence code is zero. 
Tables are input in free format as is shown on the following 
page. 


68 





4 
ct 


Contents 
Number of data pairs to be entered (This 
Should be a single card) 
Temperature 1 
Property 1 (These two entries should be on 
One cara) 
Temperature 
Sea 2 
one card) 
e etc 
N Temperature N ; 
22N+1 Propert s (These two entries should be on 
one car 
If any surfaces have a_ specified boundary condition 
(NSBC>0), a card must be submitted for each surface. If the 
same surface is subjected to more than one pommecny Conca 
tion a card must be used for each one of these conditions. 
The total number of cards must equal NSBC. The information 
is entered in free format as follows: 
Item Contents 


2 
(These two entries should be on 


#¢+ NS WH =o 


No 


1 Element nunber . 
2 Boundary condition code (see below) 
3 Surface code (see below) 
4 Property value (see below) 
5 Ambient temperature 
The ues yaar Oe codes are: 
ux 
2 Convection (constant coefficient) 
3 Radiation 
4 Convection (temperature dependent property) 
The line codes are: _. 
1 R = +1 line 
-1 R = -1 line 
2 S$ = +1 line 
-2 S$ = -1 line 
3 T = +1 line 
-3 T = -1 line 


, 
| 


oO a 





@ 


The property values are; 
Flux - Flux per unit area aa 
Convection - constant heat transfer coefficient 
or eee Lf temperature dependent) 
Radiation - Product of emissivity by Stephan- 
. Boltzman constant 
The ambient temperature is ignored for the flux boundary 
con dition. i. ; 
If the boundary condition code is 4 Geerrer eta dependent 
heat transfer eee cee cel a table mus follow for the 
temperature dependence. £ is input in free format as is 
shown on the following page. 
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Item Contents 
1 Number of data pairs to be entered (This 
should be a single card) 
2 Temperature a 
3 Heat transfer coeficient 1 (These two entries 
should be on one card) 
+ Temperature —— ; 
5 Heat transfer coefficient 2 (These two entries 
should be on one card) 
etc etc 
2 N Temperature N _ 
2eN+1 Heat transfer coeficient N (These two entries 
should be on one card) 
CTEMP 
This operation zn EUtS constant temperature boundary 
restrain ata. I ou have nodes with constant tempera- 
ture, you must enter the temperatures for those nodes. to L 
nodes with no temperature restrictions, no entries should be 
made. Automatic generaticn Capability 1s built into the 


Operation. Data is entered in free format as follows: 
Item Contents 
1 Initial node 


2 Last node 
3 Node increment 
4 Temperature 


These entries may be repeated until al} the constant temper- 
ature nodes are entered. This operation must be terminated 
by a row of alternating zeroes and blanks. 


PROF 

This operation establishes the profile of the equations for 

solution of the heat transfer .problen. After the issuance 

of this ccmmand the problem is set and you may not change 

the node mumpers with restrained poe ey conditions 
ra 


(constant EChbe panures - The values of the res ined temp- 
eratures may be changed. 
4 


This operation forms the symmetric conductance matrix for 
heat transfer problems. 


This operation forms an unsymmetric conductance matrix. 


ee soo fe oi ne 


This, operation forms a lumped capacitance approximation 
Matrix for heat transfer problems. 
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= =e apa 42 a oe ane 2P a S aoe ee ee? eae ees SFP SS PP BS 28 SS Sew SP Se SS Se eee ee eee ee ee ee ee Se eee ee ee ee eee Se ee Ss = 


CCAP 


This operation forms a consistent capacitance approximation 
matrix for heat transfer problems. 


np Een 2 ap «<p 22a ea ap Sean «42 ae 42 ea a SP 42a oP? COR ee eee OES Bees 2 ae SP SE Se Se Se Se ee Oe SEP eee ee ee eee eee Se Se ee ee Se ee ee 6 6» a2 a a» 


tes operation forms the flux vector for heat transfer prob- 
ems. 


CALC 


This operation solves time independent heat transfer prob- 
lems for temperature and updates the temperature matrix. 


Berean 42 a 22a ee 22a 2 Saw ae 224 ae a ae 22 SP SP eS ew 22a 2S SS eee 2 SE Se SS SS Se ES Se Se 2 ee Se SP SO SS SO Se ee ee ee a 2 a ee 


ODE ,M1 


This operation, solves the, first order o dinary differential 
equations arising rom time dependent heat transfer prob- 
lems. If M1i1is "INIT",the initial conditions are _ estab- 
lished. This operation must be followed by the following 
data in free format: 
Item Contents. . 
1 Integration paramter theta for two point 
scheme (default = 2/3 
2 Integration parameter gamma for three point 
scheme (default = 1.5 
3 Integration parameter beta for three point 
scheme (default = 0.8 
4 Maximum temperature difference for time step 
adjustment 
5 Ninimum temperature difference for time step 
adjustment 
ae values ate obtained by entering zeros for items 1, 2 
an . 


If M1 is "LINE", the two point eee scheme is .used. 
Tf M1 1s "QUAD", the three point integration scheme is used. 


Some suggested values for theta for the two-point scheme as 
given in Reference 4 are: 


THETA 
Crank-Nicolson 1/2 
Zlenklewicz 2/3 
Bettencourt 3/4 
Liniger .3878 


Some suggested values for beta and gamma for the three point 
scheme aS given in Reference 4 are: 


Beta Gamna 
Lees lf 3 1/2 
Hoage 3/4 2:0 
Wood . -646 1.2184 
Zienkiewl1cz u/5 3/2 
Bettencourt 9/10 3/2 
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This operation prints the nodal temperatures of a hea 
transfer system in node number order. 


= ape 427 ae 2p oP oP a ee 2 ape 2p em 2 oe om 6 22a 22 ew eee ee ee ee) ee ee ee Se Se eee ee ee se 6 2 ee 2 2 ee ee 6 2 2 2 2 2 Se 2 oe ee 2 ee 2 a oe aD 


TOL,M1 


This operation. sets the solution convergence tolerance to 
the value found in the 1x1 matrix named M1. 


CONV 

This operation performs a temperature convergence test on a 
heat transfer systen. If this operation is used inside a 
loop, (LOOP. operation) and the test shows convergence, 


looping will be terminated. 


DTIM, M1 
This operation sets the time increment for integration in.a 
heat transfer system to the value found in the 1x1 matrix 


ADT IM 

This operation advances the time in a heat transfer on ew 

by one. time step. The time step is input with the DTIM 
operation. 

EEE EEE EE 


This operation computes the dominant eigenvalue and eigen- 
vector of the current heat transfer conductance matrix. 


PROMPT 


This operation permits the ana FOE USEE LMput to be 
Baie once without loss of other output. A second PROMPT 
wi restore the prompts unless eae print output 1s 
suppressed by the NO command. 
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5. Graphics Qperations 


CAL“-NPS has a limited _graphics capability for users 
with EPLOT-10 compatible terminals or IBM 3277 Dual Screen 
terminals. fTwo- or three-dimensional meshes may be viewed 
in the X-¥, Y-Z or X-Z planes. 

The graphics operations are initialized with the 
GRAPH operation by which the user specifies which graphics 
capable terminal is being used. The PLHX operation plots 
heat transfer meshes and the PLST operation plots structural 
meshes. The user may title the plots using the TITLE opera- 
£0 De 


GRAPH 


This operation initializes the graphics , packag : Presently 
there are two terminals users May utilize f6r *he graphic 
Operations. A yes must b2 input immediately following this 
card to an that one of he two types of terminals is 
pean ceed. Then the terminal code must be entered in free 
or mat: 
Terminal Code Terminal . 
1 PLOT-10 Compatible Terminal 


Z IBM 3277 Dual Screen 

EEE EE EE 
This operation plots .2-D and 3-D heat transfer analysis 
eoes- If the mesh is 3-D, the user may speciry the plane 
S Lews 

N1 = 1 : X-¥ plane 

N71 = 2 : Y-Z plane 

N1 = 3 3: X-Z eeeee . 
The default is the Xx-Y plane, hence N1 is not required for 


2-D meshes. 
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This operat ~on plots, 2-D and 3-D structural analysis meshes. 
M1 is the node coordinate matrix created by the NODES opera- 
tion. S2 is the element connectivity matrix loaded with 
the LOADI operation. The row. dimension is fhe number of 
elements and the column dimension is the maxiaum number of 
nodes er element plus one. M2 contains the following 
informatior: 


K,1) = Number of nodes in element number K 
M2(K,2) = Node 1 
M2(K,3) = Node 2 
atc. etc. 


aera 1), = Node N . 
N14 moee? = <5 he plane of view: 


1 = X-Y plane ’ 
N1 = 2 Y¥-Z plane 
N1 = 3 #: X-Z plane . 
If the structure contains two dimensional membrane elements 
(PLANE operation), the connectivity must follow the conven- 
#zo0n shown below: S 
7 Zz 





Ny Noi8,-1) (1) 





TITLE ,N1 


This operation allows the user to input N1 (up to «hree) 
fifteen character lines to label plots generated with sither 
the PLHX or PLST operations. The label will appear in the 
upper right hand area of the screen outside the plotting 
area. . This Ree o must pate preceed the Lotzing 
operation. The title may be changed by reissuing the TITLE 
command. The default value of N1 is one. 
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6. Looping Qnerations 


CAL-NPS has a five level looping ability. fhe first 
Operation is LOOP and the last operation is NEXT. 
Operations within CAL-NPS are normally executed as they are 
encountered. If the operation requires data, the data cards 
follow the operation card. In the case of looping, however, 
all operation cards are stored within the computer before 
they are executed. If operations within the loop require 
data, the data cards must be supplied in the order required 
after the last NEXT operation. If an error is encountered 
while executing ina loop, the entire matrix of loop 
commands 1s deleted and the user is given the opportunity to 
try again. Matrices that have been modified by operations 
successfully completed while in the loop remain modified. 
After all loops are executed the computer storage required 
for these operations is automatically released by the 


program. 
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LOOP, N1 


N1 ais the number of times ee is to be executed. 
Associated with each. LOOP operation, there must. be a corre- 
sponding NEXT operation which signifies the. end of the loop 
and the return of the control to the beginning of the loop. 
The following is a possible series of looping operations. 


First level loop 





executed 5 times 


LOOP, 2 Second level loop 
aoe executed twice. 
aoe (total cf 10 times) 


LOOP, 4 Second level loop 
7 executed 4 times 
(total 20 times) 





NEXT,M1 or NEXT 


The operation NEXT signifies the end of a loop. It is 
apparent which LOOP and NEXT cards are associated if there 
are an equal number of each. The operation _NEXT,M1 will 
cause the loop to terminate if the first term in the matrix 


M1 1s negative. 


SKIP,M1,N1 


This operation will cause + 
1f the first term in the ma 


c+ 


e ski of the next Ni, operations 
rix named M1 1S negative. 


ap eS! 2222? 22 OR ee SP ae ee ae O82 ee Pee ee 22 Oe ee es eee ee See Oe 22 2 ee ee SOP 6 ee ee Se ee ee 2 ee ee ee ee ee Se ee ee ee ee 2 ee ee ee 2 Se ee oe a = op 


7. User Defined Operations 


USERA and USERB 


These names are reserved for operations to be defined and 


ayaa enmec by the user. In order to program these opera- 

10nS it 1s necessary to understand the internal organiza- 

eee: oo CAL-NPS. Chapter III of Reference 2 contains the 
etails. 
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E. LARGE PROBLEMS 


CAL-NPS is designed as an educational tool. It does not 
take advantage of banding and symmetry in matrix storage, 
except in the heat transfer operations. Larger problems can 
be solved by increasing the dimension of the L array, but a 
general purpose program that makes maximum advantage of 
out-of-core storage and takes advantage of banding and 
symmetry for in-core matrix storage is probably a better 
choice. With the above disclaimer, to increase problem size 
capability, increase the dimension of the L array and change 
the value of MAX toethe new dimension size in the following: 

G=—==-F/AIN PROGRAM 
G—===—-Sil PROGRAM CAPA CITY 
COMMON MTOT,NDP, L(100000) 
MTOT = 100000 
NDP = 2 
CALL SETIME 
CALL CAL1 
STOP 
END 


With the dimension of the L array as above, the program 
currently executes in 1024K bytes for CP/CMS. The region 
necessary for execution will increase about eight times the 


increase in the L array. 
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APPENDIX B 


SAMPLE DATA FILE 


file (FILE 


hollow cylinder with circumferential heating strips problen. 


the 


for 


FTO4FOO 1) 


Sample data 


isa 


This 


03-2. 
Ne 
80 2 1 
ace 1 
28 2 1 
16 2 1 
05 2 1 
93 2 1 
O52 1 
3412 1 
O32) 1 
4212 1 
080 2 1 
So 2. 4 
eaoee % 
116 2 1 
B05) 2 1 
493 2 1 
870 2 1 
3412 1 


ae Oo Cu ak mo ssioN cee MEO IKG DaAncOne sna moe 
mwHoowonef Dn F FY NN FTMNIMWNOMOWOMOTORnWATO0Bd\O wo 
DONSFONOTONSSON OMS Oo OO N oF N 
e§ ©C 0 oO OH #eO #O #O 2O #0 eae © 

AD FHAKDAIDANAANANKAEKE AKO rn OM FmOMrNnN NN NAN &f- & © 
en ee ee Se ee GS AO OO. Gk) OO 
Or OF OP OF OFF OP OF OF OVW COODODODODOVDOWDODOWDODVODBDOWDO 
WMO MOM HO MID MYO MOM IM OMO MOMNANANHANDHN HANAHAN Muni or 


9 
q 
Z 
0 
7 
3 
7 


WOFOKMOKOTOUKOKOKTOKTUOFrOryOKraMsgMaAnMsamMsaMsaMmnsanmnanmsansnsm 

- WOOMVOOOWOODOWOCOOCVWUOCOOVUOCVUOSOrOrOrorerOr er OerOrOroeoO 
WOWODOOODOWODODOOWUOOOWODODVUOWOOWOWODWOWOWODOWODOWOWOWO 
MOTOAWO eO eO oO eO eO etO eO e\O ee Oh Om Ohm Oh OF 8h Oh Om Om OT Om CO CC 
Ri Mer tiem Ie tiem rime im tnm di et ot oh oF of oF oF OF oF oF oF oF of 

© e a e e e e e e e 

Ho) oS -Fo- = ee F fF Ff © FOC CC TOMNTCOSE EP CDE NrorerenrMnrsae 
CDNA QS ONE MS FONT OSPR COCKS KT CSC KM KH MNOEM KNEE Mersae ses eseasoH 


78 





© 
© 
| © 
© © 
© © 
CF POS Tl © © 
NNN NNN NN NN © © 
MO ror wmdvown mm Oo & © © 
ON DOWN N wf SO HH ~~ = - FTF Tl lhe 
PMRNTOOSKMAAIMK ASMNAFODrMD MoO Oo a 
ecg SON OV De ON ae | e Te eo 
mowoodoewon wf WD st @® = © © 
DONFONSTONSOOMOMORONS ONOMOATOMOOSHrOWO 
OO TPeEANAMMATNNMOO’H DO © © 
HN SMMANIAAINANANe NT NoOMNon © 
KP FP FY FY FP FP PY FP KF PFO KFTOoOMoeroncvucronmoncoo © Oo oooocw 
MWD INWUIWWM OMUDMNOLIWMNWINWOLNO - NAN OM FF MW woO ™ © mM worovwow 
MOMOMOMNOMOMOMOMOMOAMOMaEo AY Ff LF KF KF KF FT KO DBDWO © 
EO CO KOCK TE KCK eR rT wmnwoiPnr Dn oO LO wow CDo0oo-. 
Pee rich eh eRe CR eM eR el KO MOEN KM HTM CM KEH KNEKO wee O COecrren 
zszov7vo7o70o7 0790 3070202702370 '" N Ww) 
om et oem OE Oe AO MO CP TPFEM CUCM COWDECOCONAURKHoO PLN OIOANCN ly & 
ef ef ef ef oF oF oF oF} oF eH ok ON N OM FNM OO & DW =a oO A e Nm&ORURA 
OUOere FF Ff & - & “- er OMo oToNn COIN CIC OF WMA «x 
COST HNTOCCOEN CNC TOENCTOrr&— -1 cr FF TO KE Tl ET hue H om = OMMmO at 4) 


OPPO WO CM OE HO OT HOCK NEM EK TON COFFE KOUTDOAeE @eTOrONDraAHNmUAM 
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APPENDIX C 


EXAMPLE TERMINAL SESSION 


The following terminal session was recorded using the 


data file in Appendix B. 


CAL 
ENTER TERMINAL CODE: 

1 = PLOT-1Q0 Compatible Terminal adi APHICS) 

2 = IBM 3277 DUAL SCREEN (GRAPHICS 

3 = Any Alpha Numeric Terminal (No? GRAPHICS) 
IBM 3278 TERMINAL 
CP TERMINAL LINESIZE 132 
GLOBAL TXTLIB CMSLIE FORTMOD2 MOD2EEH NONIMSL IMSLSP 
Pep eoer OF DISK FILE OT A VS LRECL 3408 BLKSIZE 3412 
FILEDEF 908 DISK FILE FTO8FO01 (RECFM FBA LRECL 132 BLKSIZ 
FILEDEF O04 DISK 
PILEDEF 13 DISK CAL TEST (RECFM FBA LRECL 132 BLKSIZE 132) 


FILEDEF 50 DISK HLP CAL C 

LOAD CALO CAL FRTCMX GROUP1 GROUP2 GROUP3 GROUPY GR (CLEAR 
EXECUTION BEGINS... 

wan == == --- = -- == - -- - -- --- --- --- 0.0 SECONDS 


START 
LAS TART 
— _ 0.003 SECONDS 
PROMPT 
a DROMPT 
epee 0.003 SECONDS 
READ, 4 
wAR BEAD, 4 
aH TX FR 
99 ROWS 1 COLUMNS 
1 ROWS 8 COLUMNS 
2 ROWS 4 COLUMNS 
4 ROWS 1 COLUMNS 
1 ROWS 4 COLUMNS 
4 ROWS 1 COLUMNS 
4 ROWS 4u COLUMNS 
7 ROWS 1 COLUMNS 
10 ROWS 1 COLUMNS 
1 ROWS 99 COLUMNS 
2 ROWS 99 COLUMNS 
5 ROWS 80 COLUMNS 
1 ROWS 99 COLUMNS 
99 ROWS 1 COLUMNS 
99 ROWS 1 COLUMNS 
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NDNOorrtd=s 
NW 


0.003 SECONDS 


NODE 


MBER OF NODAL POINTS 


NU 
NUMBER OF ELEMENTS 
DIMENSION OF COORDINATE SPACE 


/ 
NODES PER ELEMENT (MAXIMUM) 


NUMBER OF MATERIAL SETS 
# COORD 


DEGREES OF FREEDOM 


Ht NOOMMmOMIAN MNOOPFWOWONIN MNOCOFoOowOomMmaran ONOCOrFOUMAN ONOOrNoOOMAN worn 
TT OOOCOCOooO KFT rTOoCooCceeo FeKrOCOOCCOC0CO FKrCCCCo0o0o0 &FrrCOCOCCO0O xr-K-o 
DOO OOOO 9NOOOOOHODOOOOONOOOO COSCO DOCS COCO COCO CCDDCOCOCCODCCOOCCOSD 

*ee0eveseteeeeeeseseteeveeeeere eevee ewvvreeer ewe we ee we wee eh el hl hl OC hl hl OC CCl Ml hl hl tl le 


DOOD OOO OOOO OOOO OOSCO COCO COSC OHOVOCOCSCOCO COCO CODCSCCOCCCOCOCCOCCODNDO 


RODDOUEEE EEE NOUNMNMNMNNMMMNMNMN Satta eae tg Qrr rr KK PCT eK KE PEPE EEE er ODMOMMM 
(0.010 1010109119 101019191010 10110 1010 10.10 191.9110 10 10 19. 110 10.110 19 10. 10. 1010110 10.119 (L011 (10.10. 10. (0 LOLI LALA EAL LAA LAE LALA 


0 0 CF C07 CT C07 FCC CTCF CC CC CO CC CTC CC CC CFP CC CC eC eC CCC er cece 
ee | 


COO OOOO OOOO OO OOOO OO OOOO COSTS DOOCODSCCCOCCCODOCOCCOCOCOCCOCCOCCOCCOO 


TFHOANM FMF DNOSC NM PMNOF DHOK NUM PNOPDHOK AMANO DHOKNMANOFODHOKANMaANOr co 
ETS TET ET NNNNNNNNANMMMMMMOMMNMN sata aessetes II NNMMMM MM 


© 
Zz 
Fed 
A 
Oo 
= 
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-9083 SECONDS 


DBDM™OWMIN MEK A™-OWOMUIMN FOOFOMWMAIMN MANHDOOMMAIAMN 
COOOOOOD FrTrOoQo0QO009SoO FF KF DOOD OO0ooO FOQOooo000OoOoOo oOo 
welelelolelolelelelelololelelelelelelololelolelololelelelelolololejelolelololole 
OOOO ODO ODO OODOD ODODODO VOODOO ODO OO DVO OOOOOMODOVOOOO0O0O 


MSA STITATIOCL LESCOL SE LNHOOWOUWNMOUWONUNUU FTNNMMNNNMNMNMMNW 
WWW WIA) FF SS TT IMM MMM MNMMMMMANNANAANANINANAN 
Sete Geen Comte OOo Oe Oe FO Ce CS SSCS CS CCC CCC SC CSCS SK 
eee eee ere eee ee ee ewerwewe ee em ee em em hm em he OC hel Oh Ole eh 8 


eeolelelealelalelelelelelelelelcelolelaleleleleolelolololelelololololololelelololole) 


DOCNMNSTMNOPDHNOKAUMANOEDHOCAM THOM DHVLEAMANMOP WO 
MWODDODOOODODE EE EEE E EEE DDD DD DODDDAANAHAHAAHAAAH 


2RBE LCON 


TFONMAATMNOSDAHAOAM APMOP DHOKEMAMN 
Se EE EK | NIONANN 


ANAM ANNO DHNOT MANOS DHSKNANO 
SOE EH KH ONANNANNANNN 


NOAFNVOEDHDENANOl DHOENMMNMor 
SE MN AN AN ANCA MMAMMONAMMNM 


FHAMANOl DQDHOTMANOlSOHOSENANHNO 
EEE KKM NNNANNNAANANNANMMMOAMM MH 


THNMNANOEDBHNDK NAM ANOP ODHOKNM 
EE EM ENNNN 
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0.037 SECONDS 


Ol DHNOSCNAUSTMNOMP DHOSKNAUMNMNOlP DNOTAM FOP DHNOEAUMNSFWMM DHNOE AUN ANODHODE-AM 


ZuNnNWOr~ 
AQNANMMMAMMNMMM MN st SIA AAT TA AMMMNMNNMNMNMI NY OH OOWODW DOUOFE EEE ES OODDODODODM 


POND KAMWDWE DHDNOK AM FOP DHNOKAUM FONE DHNOSCAMANODNHD KE AMANOEVHDOE AM AMNWl © 
AQANAMMOAMNM OM OM) =F FTF STA TNMNNNMNNMN NWI OW OOOO DODOE EEE EEE DODDDOODO 


DHOKAUM FO DNHNOSCAM SNE DHOKAM AN ODNHO KH NM ANOK HOST ANMANOl DOK ANMAMNOM Dor 
NO SF AAAS STS MNMMNMNMN NNW OW OW OO 00 El EE EEE ODO DOD DDAOANAAHANAAHAAAHAAA 


MP DHDEKHAUMN OF DHNOTNAM FOlé°r>dDHNOCNAM TSU DHDNOE AM FNO DHE AVMANOFEAOlCAMANOM © 
MOM FSS FAS ATANMNNMNMNMNI NII 0 OW W010 OW DOE PEEP OW OD OD DOD DOHAAAAHAAHAHAHO 


FINO DHOT AM FINO DHNOSCAM FNO DHDNOECAUM ANNO DNODEK NUM ANO DHOS AM SOP DHS 
NQNANANANMMAMNMMNMN NM NM Stee PFAASATATMNNMMNMNWMNIMNMNMNWN46 OOW OW WOWOVDOEE EEE EEE hho 


22C TEMP 
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NODE TEMPERATURE 
89 60.000 
90 60.000 
97 60.000 
92 60.000 
93 60.000 
94 60.000 
95 60.000 
96 60.000 
97 60.000 
98 60.000 
99 60.000 
~-= =-------+-----+---------- --- 0.007 
PROP 
MATERIAL SET 1 FOR ELEMENT TYPE 2 
DEGREE OF FREEDOM ASSIGNMENTS 


HEAT CONDUCTION ELEMENT 


CONDUCTIVITY KX = —.1000000 
SPECIFIC HEAT -6900000 
DENS ITY 70. 00000 
HEAT GENER/UNIT VOL .0 
4 GAUSS PTS/DIR 
5 LINES WITH SPECIFIED BOUNDARY 
PLANE ANALYSIS 
5 ROWS 3 COLUMNS 
3 ROWS 9 COLUMNS 
5 ROWS 2 COLUMNS 
LINE B.C. 
ELEM B.C. LINE PROPERTY VALUE 
5 2 2 10-00000 
7 2 2 10. 00000 
8 2 2 10. 00000 
9 2 2 10. 00000 
10 2 2 10. 00000 
~-- ------------=----------- --- 0.027 
cP RO F 
82 ROWS 1 COLUMNS 
~-- =------- ---------- -+---- --- 0.017 
ams YMC 
945 ROWS 1 COLUMNS 
~=- +----=---- ---=-=-- -+---- --- 0.406 
POR M 
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